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We outline a formalism to carry out TDHF calculations of fusion cross sections for spherical + 
deformed nuclei. The procedure incorporates the dynamic alignment of the deformed nucleus into 
the calculation of the fusion cross section. The alignment results from multiple E2/E4 Coulomb 
excitation of the ground state rotational band. Implications for TDHF fusion calculations are 
discussed. TDHF calculations are done in an unrestricted three-dimensional geometry using modern 
Skyrme force parametrizations. 

PACS numbers: 21.60.-n,21.60.Jz 
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I. INTRODUCTION 

Heavy-ion fusion reactions are a sensitive probe of the 
size and structure of atomic nuclei 0]. Recent exper- 
iments with radioactive ion beams allow for the study 
of heavy-ion fusion with exotic nuclei. For instance, en- 
hanced fusion-evaporation cross sections have been ob- 
served with neutron-rich ^^^Sn beams on ^^Ni [2|- The 
synthesis of superheavy nuclei in hot and cold fusion re- 
actions 0j S 0) S 13 represents another experimental 
frontier. 

In general, the fusion cross sections depend on the in- 
teraction potential and form factors in the vicinity of the 
Coulomb barrier. Furthermore, experiments on subbar- 
rier fusion have demonstrated a strong dependence of the 
total fusion cross section on nuclear deformation 8] . The 
dependence on nuclear orientation has received particu- 
lar attention for the formation of heavy and superheavy 
elements J9| and various entrance channel models have 
been developed to predict its role in enhancing or dimin- 
ishing the probability for fusion [l^, uM ■ While this may 
be true for heavy systems, the orientation effects should 
influence the fusion process for light nuclei as well. 

There are several theoretical methods for calculat- 
ing heavy-ion fusion cross sections: a) barrier penetra- 
tion models 112. Iia. [l4l | . b) coupled-channels calcula- 
tions [ij, [13, [l3i H^; and c) microscopic many-body 
approaches such as the Time-Dependent Hartree-Fock 
(TDHF) method [11 El HI HI . 

In recent coupled-channels calculations of heavy-ion 
fusion reactions, one typ ically uses the "rotating frame 
approximation" [13, 13 which assumes that the or- 
bital angular momentum L of relative motion is con- 
served; this approximation avoids the full angular mo- 
mentum coupling and thus reduces the number of cou- 
pled channels considerably. One tends to use empir- 
ical interaction potentials which are mostly real and 
energy-independent (Woods-Saxon potential, proximity- 
type potentials, double- folding potential). An imaginary 
potential is unnecessary because one takes into account 
explicitly all channel couplings that affect fusion. Cou- 
pling potentials are usually obtained utilizing macro- 
scopic nuclear structure models (e.g. rotational model 
or harmonic vibrator). The coupled-channel equations 



are solved numerically with standard scattering bound- 
ary conditions at r — > 00, and an incoming- wave bound- 
ary condition at some point r — Rp inside the Coulomb 
barrier. The fusion cross section is obtained from the 
incoming flux at Rp- 

Fusion in the TDHF collision process is achieved when 
the relative kinetic energy in the entrance channel is en- 
tirely converted into internal excitations of a single well 
defined compound nucleus. In the TDHF theory the 
dissipation of relative kinetic energy into internal exci- 
tations is due to the collisions of the nucleons with the 
"walls" of the self-consistent mean- field potential. TDHF 
studies demonstrate that the randomization of the single- 
particle motion occurs through repeated exchange of nu- 
cleons from one nucleus into the other. Consequently, 
the equilibration of excitations is very slow and it is sen- 
sitive to the details of the evolution of the shape of the 
composite system. If one of the nuclei is deformed, one 
needs to carry out TDHF fusion calculations on a 3-D 
lattice for different relative orientations of the nuclei as 
they approach each other on a Rutherford trajectory, and 
afterwords one has to perform a suitable average over 
all orientations. Depending on the incident energy and 
impact parameter, some relative orientations may con- 
tribute to fusion while others may not. In lowest-order 
approximation, one might assume that all relative ori- 
entations of the intrinsic axis system occur with equal 
probability. This is indeed a reasonable approximation 
for relatively light nuclei where Coulomb excitation is 
negligible. However, in heavier systems, the dynamical 
alignment of the deformed nucleus as a result of multi- 
ple E2 and E4 Coulomb excitation of the ground state 
rotational band may no longer be ignored and results in 
a definite preferential alignment of the deformed nucleus 
which must be taken into account when calculating the 
fusion cross section. 

There are a very limited number of TDHF collision 
studies involving deformed nuclei, and almost all of these 
were done in an axial geometry using the "rotating frame 
approximation" , and thus could not address the orien- 
tation dependence of reaction cross sections. Recently, 
we have developed a new TDHF code which assumes no 
symmetries regarding both the collision geometry and the 
Skyrme effective interaction [23J and utilizes the modern 



Skyrme force parametrizations, including the time-odd 
terms. The numerical calculations are performed on a 
large 3-D Cartesian lattice using the Basis Spline collo- 
cation method '21'| for increased accuracy. In this work 
we shall use this code to study fusion cross sections in 
the collisions of sphcrical+deformed nuclei. 

The theoretical formalism for calculating the fusion 
cross section is outlined in Section |n] In Section IIIII 
we present fusion calculations carried out using this ap- 
proach. In Section IIVI a summary and outlook is pro- 
vided. 



II. TDHF WITH ALIGNMENT 

The evaluation of the heavy-ion collision dynamics can 
be divided into two separate steps: a) A dynamical 
Coulomb alignment calculation to determine the prob- 
ability that a given nuclear orientation occurs at the dis- 
tance r{to), where the TDHF run is initialized. The dis- 
tance r{to) is chosen such that the nuclei only interact 
via the Coulomb interaction, b) A TDHF calculation, 
starting at this finite internuclear distance r(io), for a 
fixed initial orientation of the deformed nucleus. 



bers as a function of time. In the Coulomb excitation 
code, we utilize the measured energy levels of the ground 
state rotational band in ^^^Dy up to the 18+ level |2J|. 
We also use experimental values for the reduced transi- 
tion probabihties 24] B{E2, 0+ -^ 2+) = S.SSe^fe^ and 
i3(£'4,0+ —f 4+) = 0.07e^6^ from which one can calcu- 
late all remaining E2 and E4 matrix elements within the 
collective rotor model. In this semiclassical calculation, 
one can see how the nucleus "climbs up the rotational 
band" during the collision. Negative times in Fig. [21 cor- 
respond to the incoming branch of the Rutherford tra- 
jectory, and positive times correspond to the outgoing 
branch. The distance of closest approach is reached at 
t = 0. (Alternatively, one can plot the excitation prob- 
ability as a function of the internuclear distance). The 
Coulomb excitation amplitudes, ajM(t), give rise to a 
preferential orientation (alignment). The dynamic align- 
ment formalism presented in the Appendix allows us to 
follow the nuclear alignment as a function of the internu- 
clear distance vector r(t). We would like to stress that 
the quantity which enters in our TDHF fusion calcula- 
tions is not the Coulomb alignment after the reaction 
has taken place (at t — > +cxd) but rather the alignment 
at the finite internuclear distance r{to) ~ 15fm on the 
incoming branch of the Rutherford trajectory. In the 



Coulomb excitation 



For a given incident energy Ecm and impact param- 
eter b, we carry out a semiclassical Coulomb excitation 
calculation of the dominant ground state rotational band 
of the deformed nucleus (see Fig. ^). The Coulomb ex- 
citation calculation starts at very large internuclear dis- 
tances (about 1500 fm) when both nuclei may be pre- 
sumed to be in their respective ground states. We have 
provided the details of this formalism in the Appendix. 
Fig. [21 shows a representative result in medium-mass nu- 
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FIG. 1: Coulomb excitation of the ground state rotational 
band members via multiple E2 and E4 transitions. 
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'^62 ^Vj where we plot the excitation proba- 
bility of the ^^^Dy ground state rotational band mem- 
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FIG. 2: Coulomb excitation probability of the 0"*", ..., S"*" 
members of the ground state rotational band in ^^^Dy as a 
function of time. The deformed nucleus is excited with a ^"^Ni 
beam at Ecm ~ 200.1 MeV at impact parameter b — 0. 



case of a spherical-|-deformed system, the relative orien- 
tation of the deformed nucleus may be specified by the 
three Euler angles (a, /3, 7) relative to the collision plane 
as shown in Fig. [3 We use here the definition of the 
Euler angles given in Ref. |2y|. Multiple Coulomb exci- 
tation predominantly excites the members of the K — 
ground state rotational band which preserves an axially 
symmetric shape; hence, rotations about the symmetry 
axis z' described by the Euler angle 7 are irrelevant (see 
Appendix). 




FIG. 3: The orientation of the intrinsic body-fixed frame of 
the deformed target nucleus is specified by the three Euler 
angles {a, l3, 7) relative to the collision plane. 



B. Fusion cross section 

Subsequent to the determination of the initial condi- 
tions for the two nuclei, as described in the previous sec- 
tion, we perform a TDHF calculation to establish the 
outcome of the collision process. In TDHF, for a given 
set of initial conditions (energy, impact parameter, orien- 
tation, etc.), only one outcome is possible, thus we define 
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where b is the asymptotic impact parameter, Ecm is the 
associated center of mass energy, and quantities /3 and a 
denote the orientation angles discussed above. The total 
fusion cross section may be written as 
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The probability for fusion is given in terms of the differ 
ential probabihty 

^-^ fusion K'^-) -^cm-) Pi ^j 
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where the Euler angle solid angle element is given by 
dn = sin(3d(3da, and 

dP fusion {b, Ecm', P, Q;) _ dPorientjb, i?cm,^0;/3, Ct) 

dh ^ dh 

X PTDHFib,Ecm;0,a) . (4) 

Since for a given Ecm value and angles (/3, a) the TDHF 
fusion probability is either zero or unity we can instead 
write the integrals over the impact parameter to termi- 
nate at bmaxiPjo), where bmax{l3,a) denotes the maxi- 
mum impact parameter for fusion for orientation angles 
/3 and a. This allows us to write 
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sin^dpda 

The differential cross section given in Eq. jSJ) in terms of 
Euler angles should not be confused with the laboratory 
differential cross section given in terms of the scattering 
angles. The differential orientation probability used in 
Eq. |SJ| is evaluated in the Appendix, Eq. (|A.7p . 



III. NUMERICAL STUDIES 



In order to gain a better insight into the differential 
orientation probability mentioned above we have per- 
formed calculations for systems that involve two heavy 
reaction partners and for systems which are composed of 
two light nuclei. We first examine the dynamical orienta- 
tion, due to multiple E2/E4 Coulomb excitations, of the 
deformed nucleus ^fgDy colliding with a spherical IsNi 
nucleus at Ecm — 265 MeV. In a central collision (im- 
pact parameter 6 = 0), the orientation probability of an 
even-even nucleus is independent of the Euler angle a 
because, starting from the 0+ ground state, only mag- 
netic substates M = of the ground state rotational 
band are excited. Therefore, the orientation probabil- 
ity depends only on the Euler angle (3 in this case (see 
Appendix, Eg. IA.7|) . Fig. 0] shows the orientation prob- 
ability for two values of the internuclear distance. At 
very large distance {R = 1363 fm) before the collision 
where the deformed nucleus is in its 0+ ground state, all 
Euler angles (3 occur with equal probability, and there 
is no preferential alignment. However, at an internu- 
clear distance oi R = Ri + R2 + 2 — 13.2 fm (where 
i?i and i?2 denote the mean square nuclear radii), we 
observe substantial alignment preference. Clearly, the 
perpendicular orientation (/? = 90°) is preferred over the 
parallel (/3 = O'^) orientation, with an alignment ratio 
dP{(3 = 90°)/dn I dP{P = 0°)/dn = 1.54. The reason 
for this preferred alignment is easily understood in terms 
of classical electrostatics: for a given internuclear dis- 
tance i?, the perpendicular orientation of the deformed 
nucleus minimizes the Coulomb interaction energy which 
is proportional to P2{cosj3). 
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FIG. 4: Dynamic alignment due to Coulomb excitation of 
^^^Dy by a ®*Ni beam. Shown is the orientation probability 
as a function of the Euler angle /9 in a central collision at 
internuclear distances R = 1363 fm (blue curve) and at _R = 
13.2 fm (red curve). 



We have repeated the above analysis to study the 
Coulomb excitation of the deformed nucleus IgNe in a 
central collision with an ^|0 nucleus. In the Coulomb 
excitation code, we utilize the measured energy levels of 
the ^^Ne rotational band up to the 6+(6.311Mey) level 
m. The B{E2,0+ -^ 2+) value is computed from the 
measured halflife J25j | T1/2 = 3.63ps, and the remaining 
E2 matrix elements are determined from the collective 
rotor model. Fig. |S1 shows the corresponding orienta- 
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FIG. 5: Dynamic alignment due to Coulomb excitation of 
■^■^Ne by a ^^O beam. Shown is the orientation probability 
as a function of the Euler angle /3 in a central collision at 
internuclear distances R = 154 fm (blue curve) and at i? = 14 
fm (red curve). 



tion probability for these two light nuclei. Despite the 
much higher incident energy oi Ecm = 95 MeV (about 6.9 
times the Coulomb barrier height), we observe only a very 
small alignment preference: dP{f3 — 90'^)/dfl / dP{f] — 
0°)/dri = 1.007. There are two aspects pertaining to the 
importance of dynamical alignment in calculating fusion 
cross sections. In the case where all orientations lead 
to fusion the alignment probability provides a properly 
weighted sum for calculating the fusion cross section us- 
ing Eq. ^. However, implicit in the above expression is 
also the fact that for different orientations we will have 
different values for bmax- Thus, even in a case where the 
central collision leads to fusion for all orientations of the 
deformed nucleus, and even if the Coulomb excitation 
probabilities are angle- independent, we still have differ- 
ent contributions to the fusion cross section arising from 
different orientations. 

In order to demonstrate the points made above we have 
performed TDHF calculations for the ^^0+'^'^ Ne system 
at Ecm = 95 MeV. We have chosen this light system at 
a relatively high beam energy so that we can perform 
the numerical calculations faster. We have used our new 
TDHF code [23, which works in 3-D and uses modern 
Skyrme forces, including the time-odd terms. For this 
work we have used the SLy4 parametrization [27|. In 



this case Hartree-Fock calculations result in a spheri- 
cal ^^O nucleus and a ^^TVe nucleus with strong axial 
deformation. When no spatial symmetries are assumed 
Hartree-Fock calculations generate an orientation for the 
intrinsic coordinate system, with respect to the code co- 
ordinate system, typically determined by the choice for 
the initial single particle states, e.g. harmonic oscilla- 
tors. Depending on the order in which the Cartesian 
oscillator shells are filled we get a particular orientation 
for the nucleus. Different orientations with Euler an- 
gle (3^0 can be generated by rotating the coordinate 
frame in which the initial single particle states are cre- 
ated with respect to the code frame. The resulting states 
will then be oriented with Euler angle /3 with respect to 
the code frame. A subsequent rotation perpendicular to 
this direction would then generate the a rotation. The 
Hartree-Fock iterations in generating the static solutions 
preserve this orientation since in an unrestricted geome- 
try all orientations are exactly equivalent. 

We have performed calculations by keeping the angle 
a fixed at a = 0° and varying the angle /3 in 10° inter- 
vals. At this energy, all angles up to f3 — 60° do not 
show any fusion for head-on collisions, whereas larger 
angles fuse. For each of the orientations leading to fu- 
sion we have made a sweep over the impact parameter 
to find the maximum impact parameter for fusion. This 
was initially done in 1.0 fm steps and reduced down to 
0.05 fm when the approximate fusion boundary was lo- 
cated. For large impact parameters the system undergoes 
multiple revolutions before fusion or deep-inelastic colli- 
sion. Starting from an initial separation of 15 fm we have 
followed most of the collisions to about 1200 fm/c. An- 
gles smaller than /3 = 60° contribute to the deep-inelastic 
channel. For these orientations a study of the central 
collisions show a gradual decrease in the translational 
kinetic energy between the two final fragments, ranging 
from almost 30 MeV for /3 = 0° to about 11 MeV for 
P ~ 50°. We also observe particle transfer between the 
fragments, in most of the cases the final fragments seem 
to have exchanged one neutron from the Neon to the 
Oxygen. In the fusion regime we have found maximum 
impact parameters for fusion as 6.35 fm, 6.55 fm, 6.83 fm, 
and 6.87 fm for f3 values of 60°, 70°, 80°, and 90°, re- 
spectively. In Fig. El we show various time slices of the 
TDHF evolution for the case of /3 = 60° and b = 6.35 fm. 
During this time interval the system makes about four 
revolutions and eventually the internal structure shows 
no memory of the initial structure. 

For Euler angles (3 = 60°, 70°, 80°, and 90°, re- 
spectively, we find ^^O -f'^^ A^e fusion cross sections 
d<J fusion{.Ecm = 'd'oMeV]l3,a — 0°)/ sin(3dj3da of 
633 mb/sr, 673 mb/sr, 732 mb/sr, and 741 mb/sr. 



IV. 



CONCLUSIONS 



To date most TDHF calculations have been limited to 
the study of collisions involving spherical systems. Fur- 
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FIG. 6: TDHF time-evolution for the ^'^Ne -\-^^ O collision at an impact parameter of 6 = 6.35 fm and initial Neon orientation 
angle /3 = 60° using the SLy4 interaction. The initial energy is Ecm ~ 95 MeV. During the evolution the combined system 
makes four revolutions. 



thermore, approximations made about the collision ge- 
ometry (reaction-plane symmetry etc.), in order to make 
numerical computations tractable, have lead to the ex- 
clusion of deformed systems from such studies. In this 
paper, we have presented TDHF calculations involving 
spherical-|-deformed systems without any assumptions 
regarding collision geometry and using the full form of 
the Skyrme interaction. The details of our new TDHF 
code can be found in Ref. J3]. Dealing with deformed 
nuclei necessitates an approach to determine the orien- 
tation of the deformed system prior to the start of the 
TDHF calculations, as such calculations are initialized at 
relatively small separations. This alignment is inherently 
related to multiple E2/E4 excitations due to the Coulomb 
interaction between the two nuclei. We have outlined an 
approach for calculating this dynamical alignment prob- 
ability and have shown how to incorporate it into the 
cross section calculations. While the alignment proba- 
bilities for various orientations do not vary substantially 
for light systems they show considerable preference for 
particular orientations of heavier systems. On the other 
hand, we have also shown that the alignment has a ma- 
jor consequence for TDHF calculations. We have ob- 
served that in the ^^TVe -V^^ O system alignments close to 
the collision axis (/3 < 60°) results in no fusion whereas 
perpendicular alignments lead to fusion. Furthermore, 
the alignment naturally affects the impact parameter de- 
pendence of fusion for different orientations. With the 
advent of computer technology and numerical methods 



unrestricted TDHF calculations are becoming more and 
more feasible. It is our goal to pursue and improve such 
calculations in order to provide a tool for studying heavy- 
ion collisions with as few computational restrictions as 
possible. 
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APPENDIX: DYNAMIC ALIGNMENT OF 

DEFORMED NUCLEI DUE TO COULOMB 

EXCITATION 

We present a brief summary of the theory of dynamic 
alignment of deformed nuclei during a heavy-ion collision. 
The alignment results from multiple E2/E4 Coulomb ex- 
citation of the ground state rotational band. Details, in- 
cluding a discussion of strong nuclear interaction effects 
can be found in Ref. [23 . The associated formfactors for 
inelastic Coulomb excitation and strong nuclear excita- 



tion (proximity potential) have been derived for collective 
rotations and surface vibrations in Ref. [291] . 

There exists an extensive literature on the 'semiclas- 
sical' theory of multiple Coulomb excitation of heavy 
ions |3y|. In this approach, the excitation process is de- 
scribed quantum mechanically, and the relative motion 
of the nuclei is treated by classical mechanics. The to- 
tal Hamiltonian consists of the free Hamiltonian of the 
target nucleus, Ho{X), and of the couphng potential for 
inelastic Coulomb excitation, Vc{X,r{t)). The latter de- 
pends on the intrinsic coordinates of the target (X) and 
on the classical relative trajectory r{t). 

The Coulomb excitation process is determined by the 
time-dependent Schrodinger equation 

[HoiX) + Vc{X, r(i))] ^(X, t) = *ri^^(X, t) . (A.l) 
The eigenstates of the target nucleus are given by 



Ho{X)MX) = ErMX) 



(A.2) 



We expand the time-dependent wavefunction in terms of 
the stationary eigenstates (f>r of the unperturbed Hamil- 
tonian Hq 



^(X,t)=^a,(i)0,(X)e- 
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This system of differential equations is solved numeri- 
cally by a combination of the fourth-order Runge-Kutta 
method and the Adams-Bashforth-Moulton (predictor- 
corrector) algorithm [Sll. For the classical relative mo- 
tion r(t) we utilize the "coordinate system B" defined in 
Ref. in), page 47. 

In the specific application of this formalism to dynamic 
nuclear alignment, we describe the free Hamiltonian and 
the corresponding wavefunctions in terms of the collec- 
tive rotor model. The degrees of freedom are the three 
Euler angles X — (a, /3, 7) 



H^{X)=Trot{X) . 



(A.5) 



In deformed even-even nuclei, the ground state rotational 
band has an intrinsic total angular momentum projection 
K = Q] therefore, the collective rotor wavefunctions |2a| 
are independent of the Euler angle 7 which describes a 
rotation about the intrinsic symmetry axis z' (see Fig. ^ 



C^r{X) 



(2^)-^/2rjM (/?,«) . 
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resulting in a system of linear coupled differential equa- 
tions for the excitation amplitudes as a function of time 

ihar{t) - ^a,(i)((^,(X)|Vb(X,r(t))|0,(X)) 



The probability density at time t to find the nucleus ori- 
ented with given Euler angles X — {a, (3) is given by 
\il^{X,t)\'^; by integration over 7 we find the correspond- 
ing differential orientation probability 
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